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ABSTRACT 

Stellar-mass compact binaries in eccentric orbits are almost guaranteed sources of 
gravitational waves for Laser Interferometer Space Antenna. We present a prescription 
to compute accurate and efficient gravitational-wave polarizations associated with 
bound compact binaries of arbitrary eccentricity and mass ratio moving in slowly 
precessing orbits. We compare our approach with those existing in the literature and 
present its advantages. 
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1 INTRODUCTION 

It is expected that the Laser Interferometer Space Antenna 
(LISA) will usher in a new era for gravitational-wave (GW) 
astronomy. The galactic stellar-mass compact binaries are 
highly promising sources for LISA. These are also excellent 
sources to do astrophysics with LISA, which requires ac- 
curate extraction of astrophysical informations from gravi- 
tational waves they emit; the real goal of GW astronomy. 
Accurate extraction of astrophysical quantities associated 
with compact binaries is possible, in principle, as the dy- 
namics of these sources, when gravitational waves they emit 
enter LISA's bandwidth, is very accurately described by the 
post-Newtonian (PN) approximation to general relativity. 
The PN approximation allows one to express the equations 
of motion of a compact binary as corrections to the Newto- 
nian equations of motion in powers of (v/c) 2 ~ Graj (c 2 R), 
where v, m and R are the characteristic orbital velocity, the 
total mass, and the typical orbital separation of the binary, 
respectively. 

An important feature of stellar-mass compact binaries, 
relevant for LISA, consisting of neutron stars, stellar-mass 
black holes or a mixture of both may be that they will have 
non-negligible eccentricities, as demonst rated by several as- 
trophysically mot i vated investigations jBenacauistal 120021 : 
iGusev et alJ 120021 : IChaurasia fc Bailesl 120051) . Let us first 
consider a recent investigation that clearly demonstrated 
that a natural consequence of an asymmetric kick imparted 
to neutron stars at birth is that the majority of neutron- 
star (NS) binaries should possess highly eccentric o rbits 
JChaurasia fc B ailed Eooilh IChaurasia fc Bailesl (l2005h also 
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pointed out that their conclusions are applicable to black- 
hole-neutron-star (BH-NS) and black-hole-black-hole (BH- 
BH) binaries. Further, the observed deficit of eccentric short 
period binary pulsars was attributed to selection effects, 
in pulsar surv eys, mainly due to the ac celerated decay of 
these binaries (IChaurasia fc B ailcs 2005). In another inves- 
tigation, assuming a stationary distribution for NS-NS bi- 
naries in the galaxy, it was argued that LISA will see sev- 
eral NS-NS, NS-B H and BH-BH binaries in eccentric orbits 
jGusev et al]|2002l) . In another investigation that employed 
Monte Carlo simulations to model galactic globular clus- 
ters, it was observed that LISA may see several stellar-mas s 
BH binaries in highly eccentric orbits jBenacauistal ^002). 
Therefore, it is desirable to have accurate and efficient GW 
templates for stellar-mass compact binaries in eccentric or- 
bits. This is because LISA observations of such compact bi- 
naries shou ld provide at least the total mass of these binaries 
JSetj200ll), angular resol ution sufficient to identify their lo- 
cation feenacau ista||200 2l) and , in principle, a bound for the 
mass of gravitons jjonedboOot) . 

In this paper, we provide accurate and efficient GW 
polarizations which are restricted to the quadrupolar 
order, H+\q and h x \q, associated with compact binaries, 
modeled to consist of non-spinning point masses, moving 
in precessing eccentric orbits. These templates, which 
should be useful for LISA, are accurate as we employ 
the fully 1PN accurate orbital motion. They are efficient 
as we provide efficient numerical prescriptions to obtain 
both time and frequency domain representations for the 
relevant h+ and h x - We restricted the orbital motion to be 
fully 1PN accurate and neglected higher PN contributions 
for the followin g two reasons. The fi r st on e is that the 
inves tigations, ( More no-Garrido et a l] 119951: iPierro et all 
l200l|j . that we want to improve (both in accuracy and 
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efficiency) do not (and can not) even treat the orbital 
motion in a fully 1PN accurate manner. The second reason 
is that for stellar-mass compact binaries in eccentric orbits, 
whose orbital frequency is ~ 10 -3 Hz, it is more important 
to include the frequency shift due to the periastron advance 
(which appears at the 1PN order) compared to the shift 
caused by the rad iation reaction, appearing at the 2.5PN 
order JSetd 1200 ll) . Moreover, the methods we introduce 
here to construct numerically efficient templates can be, 
in principle, incorporated into the phasing formalism that 
gave h+ and h x with fully 3.5PN accurate orbital mot ion 
jDamour et"al"ll2004t iKonigsdorffer fc Gonakumarll200rl l . 

The paper is organized in the following way. In Sec- 
tion we introduce relevant expressions for the 1PN ac- 
curate orbital dynamics and the associated h+ and hx, re- 
quired to construct the templates. Sections [3] and 2] deal 
with accurate and efficient ways of obtaining the time and 
frequency domain GW templates. The pictorial depictions 
of our results and associated detailed discussions are given 
in Section |SJ followed by a brief summary. 



where n is the 1PN accurate mean motion, defined by 
n = 2ir/P, P being the orbital period, and et is the ec- 
centricity associated with the 1PN accurate Kepler equa- 
tion (KE) displayed below, and £ stands for Gmn/c 3 . The 
angular motion, defined by <f) and (j>, is given by 



cf>(\,l)=\ + W(l), 

X= (l + fc)n(t-to)+0o : 



(3a) 
(3b) 



where to and 4>o are some initial time and associated orbital 
phase, respectively. The expressions for W and k are given 
by 



2/3 



W{1) = (v-u + e t sinu) (1 + ^ \ 

3^2/3 



k : 



=2 ' 



v = 2 arctan 



1 + e <* \ 1/2 ^ it 

- tan — 

1 - e^J 2 



(3c) 
(3d) 
(3e) 



where 



2 PN ACCURATE INPUTS FOR TEMPLATES 

The two independent GW polarizations for non-spinning 
compact binaries moving in non-circular orbits, due 
to the dominant qu adrupolar contributions, available in 
iDamour et alJ i2004h . read 

, GmqC_ UGm 



hx\Q(r,(t>,r,4>) = - 



h+\Q(r,(p,r,(j)) = - 



c 4 R' 
2rrcj> cos 2cj> 
Gmr) 



V r 



2 Y2 .2 
+ r cp — r 



sin 2cj> 



(la) 



c 4 i?' 



+ r 



x cos 2<f> + 2rr<j> sin 2<f> 



+ s- 



Gm 



(lb) 



In these expressions, the symmetric mass ratio 
r\ — m\m,2/m 2 , mi and iri2 are the individual masses with 
m = mi + mi, R' is the radial distance to the binary and 
S and C stand for cosi and sini, respectively, i being the 
orbital inclination. The dynamic variables r and <j> denote 
the relative separation and the orbital phase of the binary 
in a suitably defined center of mass frame, with r = 4? 
and (j> = [see lDamour et ail (J2fl04) for our convention]. 
In order to obtain a prescription that models the temporal 
evolutions for /i+|q and h x \Q, namely the GW phasing, we 
invoke the following parametric descriptions, involving the 
eccentric anomaly u, for the dynamical variables present in 
Eqs. Q. 



2/3 



(Gm\ L <%, . r 

r= — (l-e t cosu)n + 

\ n 2 J [ 6(1 — et cosu) 



+ 2rj — (6 — 7v) e t cosu] 



-18 
(2a) 



6( 



e t (Gmn) 1/3 r 2/3 ] . 

T r 6 + i (6 It]) sm u , (2b) 

11 — et cosu) L J 



e = et[l + C 2/3 (4-v)} ■ 
The first time derivative of cj> reads 

£2/3 



(3f) 



(1 — etcosw) 2 \ (1 — ef)(l — e* cosu) 
x [3- (4-77)ef + (l-77)e t cosu] j. (3g) 

W e note that derivatio ns of these expressions, available 
in D amour et al. require the 1PN accurate quasi- 

Kepler ian parameterization for co mpact binaries in eccentric 
orbits iDamour fc DerueUdll985h . The explicit time evolu- 
tion for H+\q and h x \q is achieved by solving the 1PN ac- 
curate KE, present in the 1PN accurate quasi-Keplerian pa- 
rameterization, which reads 



I = n (t — to) 



e t sm u , 



(4) 



where / is the mean anomaly. Note that Eq. Q is struc- 
turally identical to the classical (Newtonian accurate) KE, 
only if we express the PN accurate dynamics in terms of 
e t , one of the three eccentricities that appear in the 1PN 
accurate quasi-Keplerian parameterization. This allows us 
to adapt the most efficient and accurate ( numerical) way 
of solving the classical KE We introduce 

Mikkola's solution in the next section. 



3 MIKKOLA'S SOLUTION TO THE 
RELEVANT KEPLER EQUATION 

The celebrated KE has enticed several generations of distin- 
guished scientists roughly from 16 50 onward and, therefore, 
a plethora of solutions exists [see IColwelll Jl993f) for a de- 
tailed review] . Even though with the help of fast computers 
it is easy to obtain accurate and quick solutions to the KE, 
there exist subtle points associated with the accuracies and 
efficiencies of various numerical solutions. A numerical so- 
lution to the KE usually employs Newton's method which 
requires an initial guess uo that depends on I and et- A num- 
ber of iterations will be required to obtain an approximate 
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solution that has some desired accuracy. The number of it- 
erations to reach this accuracy naturally depends on Uo, et 
and Z. It is therefore compelling to ask, as done for the first 
time by E. Schubert in 1854, if there exists a clever choice 
for Mo, such that with one iteration of Newton's method one 
would achieve the desired accuracy for all e« and I. 1 In 1987, 
Seppo Mikkola devised a clever empirical procedure which 
gives u with a relative error ~ 10" 15 , that may be treated 
to be one of the closest solutions to the puzzle posed by 
Schubert. In this paper, we employ Mikkola's simple and 
robust solution to solve the 1PN accurate KE and, there- 
fore, obtain h+\c}(t) and h x \Q(t) with 1PN accurate orbital 
evolution. As Mikkola's solution is crucial to obtain highly 
efficient 1PN accurate phasing, let us briefly describe below 
Mikkola's procedure. 

In order to obtain an efficient prescription for uo (an 
initial guess for u) , Mikkola introduced an auxiliary variable 
s — sin(w/3), allowing one to write our KE as 



3 arcsin s — et (3 s — 4 s 3 ) = I . 



(5) 



Truncating arcsin s to the third order leads to the following 
approximate KE 



3(l-e t )s + (4e t + -) S 3 =/. 



(6) 



The solution to the above cubic equation in s can be ex- 
pressed as 



a 

s = z . 

z 



(7) 



where 



l-e t 
le t + 



i ■ 

2 



z= [p ± + < 



1/3 



(3 being 



¥ 



4e t + i 



(8b) 
(8c) 

(8d) 



and the sign of the square root is to be chosen to be that of Z. 
Mikkola negated the largest error, occurring at I = n, by a 
simple correction term ds = — 0.078 s°/(l + et). With the 
help of Eq. (|SJ, the initial guess for u, namely uq, becomes 



uo = I + e t (3u — 4ll> ) , 



(9) 



where to — s + ds — s — 0.078 s 5 /(l + e t ). It is important 
to note that uo, which enjoys maximum relative error that 
is not greater than 2 x 10 -3 , is obtained without employ- 
ing a single evaluation of trigonometric functions. Following 
Mikkola, our transcendental equation, namely the 1PN ac- 
curate KE, also requires evaluation of just one square root 
and one cubic root for its solution. To improve uo, one needs 
to app ly a fourth-order extension of Newton's method, avail- 
able in lDanbv fc Burkardtl jl983t) . This requires us to define 
f(u) = u — et sin it — Z and its first four derivatives with 
respect to u, evaluated at u = uq. This leads to a solution 



to our KE that has a relative error not greater than 10 15 , 
which reads 



U = Uq + 114 , 



(10) 



where 



Ul = 


/ 






(lla) 


U2 = 


/ 

/' + \ f" 


Ul ' 




(lib) 


U 3 = 




/ 




(11c) 


r + \ r 


u 2 + 




114 = 








(lid) 


r + \ r 


u 3 + 





It is important to note that the method requires reduc- 
tion of I into the interval — tv ^ I ^ n in order to make s as 
small as possible. We e mploy geometrical interpret ations of 
u and I, and results from lDanbv fc Burkardl) il983l) . to map 
any / into the above interval. First consider the case where 
I < 0. We then define I* = —I and solve u* — et sin u* — 1* =0 
for u* and naturally u = — u* . If I is such that n < I < 2n, 
then there exists u(l) — 2tt — u(2n — I), which follows from 
the geometrical construction of u and I. If I > 2n, the situ- 
ation is subtle and we define Z* = I — I j j 1 2tt , where || ^r|| 
denotes the integer part of Z/27T, so that Z* lies in the inter- 
val [0,2 71-]. This leads to u(l) = jj ±\\ 2tt + u*(Z*), which is 
based on the fact that at the fixed points of KE, u — I. 

Let us now point out another way of tackling the largest 
error associated with I = ir, without resorting to the in- 
troduction of ds. 2 The idea is to choose the coefficient of 
s 3 in the Taylor expansion for arcsin s, such that Eqs. 10 
and JljJ become identical at Z = n. This can be achieved 
by replacing 1/6, the coefficient of of s 3 in the Taylor ex- 
pansion, by 1/(6 + 7 1). It is then straightforward to deduce 
that 7 = -(3/4tt)(27V3 - 16tt)/(3\/3 - 2tt) at Z = n. This 
approach provides uq with a relative accuracy ~ 10~ 3 for all 
et and as expected, there is no need to introduce ds. 

It is important to realize that Mikkola's solution only 
demands the solution of a cubic polynomial that involves 
no trigonometric functions and one-time evaluation of sin a; 
and coslj. Further, it applies to all et and Z with ^ et ^ 1. 
With the help of Mikkola's solution to 1PN accurate KE, 
Eq. P). and employing Eqs. (QQ, we obtained accurate and 
efficient temporal evolution for H+\q and Jix\q. In the next 
section, we present a way of expressing h+\Q(t) and ftx |q(£) 
that easily reveals their spectral contents. 



4 FOURIER SERIES REPRESENTATION FOR 

H X \ Q AND H+\ Q 

In this section, we present a Fourier series expansion of 
h x \Q{t) and /i+|q(i), that easily allows one to compute the 
associated po wer spectrum. This sectio n is influenced by 
Section. Ill of lOopakumar fc IveJ l|2002F) [However, we cor- 
rect few shortcomings present in that section]. We begin by 



1 See page 98. IColwefll <1993t) 



2 We thank Seppo Mikkola for pointing out this (unpublished) 
esthetically and numerically better prescription for uq. 
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rewriting Eqs. Q as 

h*\ Q (t)=^-H x \ Q (t). 
h+\ Q (t) =^ H+\ Q (t) : 

where H x \q and H+\q symbolically read 



(12a) 
(12b) 



H x \ Q (t) =X 2C (l) cos2\ + X 2S (l) sin2A, (13a) 
H+\ Q (t) =P 2C (l) cos2A + P 2 s(0 sin2A + P (0- (13b) 

We note that X 2C (l), X 2 s(l), P%c{l), Pts[l) and P (l) are 
2 7r-periodic implicit functions of / (and explicit functions 
of u) and can be obtained with the help of Eqs. Q after 
using the <j> = A + W split. Their exact expressions, in terms 
of the dynamical variables, are given by 



x 20 (/) = -¥^— 



X 2S (l) 



c 2 [\ r 
2rr^>cos2W } , 



+ d> 2 r 2 - f 2 I sin 2W 



2C ( (Gm 



+ 6 2 r 2 — r 2 I cos 2W 



P2C{1) = 



c 2 [ y r 
+ 2rr</>sin2W 

(1 + C 2 ) \ (Gm 



(14a) 



(14b) 



V r 



+ S 2 r 2 - f 2 I cos 2W 



+ 2rr</>sin2W 



P.sil) = ^l{(^ 



c I \ r 

- 2rr0cos2W j , 

Po(0 = -Cl— -4> 2 



+ d) 2 r 2 - f 2 I sin 2W 



(14c) 

(14d) 
(14e) 



Recall that the dynamical variables appearing in Eqs. I|14fl 
are parametrically given, in terms of it, by Eqs. and 
and require the solution to the 1PN accurate KE. The 2 n- 
periodic functions X 2 c{l), X 2 s(l), P 2 c{l), P 2 s(l) and Po{l) 
allow Fourier series expansions that will be employed to ob- 
tain the spectral content of /i x |q and H+\q. Let us first 
consider H x \q. The fact that we can express X 2 s{l) and 
X 2C {1) as 



+ oo 



(15a) 
(15b) 



leads to 



(16) 



where 



Si = 



(17a) 



-i2c/. 



{Cj+iS}), 



(17b) 

(17c) 
(17d) 



wt = ( 7 - + 2p), u~ = {j-2p), 
p=(l + k). 

Following iGopakumar fc Iverl i2002T) . the Fourier series for 
/fx | Q, relevant for computing its power spectrum, simplifies 



to 



H X \ P Q B = V2\ ^{S^t 1 + Cje^ 1 ) + S e l2pl 



(18) 



where Hx|q S denotes an expression required to compute the 
power spectrum of H x \q(1). In order to arrive at Eq. UHJ 
from Eq. 1161 . we have used identities like \Sj\ 2 = | C7— j | 2 and 
|S_j| 2 = |Cj| 2 . These relations are valid because 4>o, X 2 q 
and X 2 s are real numbers. These inputs give us following 
expression to compute the "one-sided power spectrum" for 
H x \q 

{OO 



* li-2p| i\ 



+ S e 



i 2 pl 



(19) 

From the above expression, we see that spectral lines ap- 
pear at frequencies (1 + 2p) /,,, (2 + 2p)f r , (3 + 2p)f r ,... 
with strength ~ \Si\ 2 , IS^I 2 , |5*3| 2 , ... respectively, where 
f r = n/2-K. Similarly, the u~ part of Eq. 1191 creates spec- 
tral lines at frequencies |(1 — 2p)\f r , |(2 — 2p)\f r ,... with 
strength ~ |<?i| 2 , |C2 1 2 respectively. Further, there will 
be an additional line at 2pf r with strength ~ \Co\ 2 = \So\ 2 - 
Let us now sketch a similar analysis for H+\q and its 
Fourier series expansion reads 



H+\q(1)= J2 (S+e^P + Cfe^'+Pfe 



(20) 



where Sj and Cj are defined in a manner similar to Sj 
and Cj, involving Fourier coefficients of P 2 s and P 2 c- In ad- 
dition, Pj denote Fourier c oefficients of Pp(l). After f ollow- 
ing the details presented in IGopakumar fc IveJ (2002), we 
obtain H+\^, an expression required to compute the "one- 
sided power spectrum" for -H+|q as 



H 4 



U+2p) I 



3=1 



+PU i3l ) + Ct e' 



+ C- e 



+ PS 



(21) 



It is important to note that in the case of H+\q there exist 
spectral lines at /,., 2 f r ... that are unaffected by the advance 
of periastron parameter k. This follows from non-A terms in 
H+\q. 

It is now straightforward to compute the one-sided 
power spectrum for H x \q and H+\q with the help of 
Eqs. 1191 and I2H . However, there are few technical details 
that need to be explained and let us, for simplicity, focus on 
H x \q. It is clear from Eq. 1191 that we require to evaluate 
Cj and Sj numerically. We employ Fast Fourier Transform 
routines from Numerical Recipes to evaluate Cj and Sj af- 
ter obtaining X 2 s(l) and X 2 c(l). Recall that X 2 s and X 2 c, 
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as given by Eqs. I14all and <14bL are explicit functions of 
u and we use Mikkola's method, described in the previous 
section, to obtain X2s{l) and X2c{l) numerically. Further, 
we use the following exact relation for v — u that appears in 
the expression for W(u) 

_i/ /3<f, sinu \ , 

v — u — 2 tan — — , (22) 

\ 1 — p<$, cos u J 

\-J\-t- 2 I 1 

where P4, = — ^— — - [See IKonigsdorffer fc Gopakumarl 

(2006) for a detailed derivation of the above relation]. It 
is in the evaluation of X2s(l) and X2c(l) that we cor- 
re ct and improve sho r tcomi ngs of a similar analysis done 
in lOonakumar fc Iverl l)2002l) . We observe that the PN ac- 
curate parametric expressions for r, r and tj> appearing 
in Eqs. 114H wer e expanded into amplitud e corrections to 
h x \ q and h j. \ n in | G opakumar fc Iyei] J2OO2T1 [see Eqs. (2.23) 
- (2.27) in iGopakumar fc Iveil (12002 ^1 . In this paper, we 
avoid this approach and treat r, r and <j> i n the same 
footing like W(u) as done in iDamo ur et all J2004T ) and 
IKoni gsdorffer fc Gopakumarl i2006T ). fri TGopakimiajfc Iveil 
( 2002), an inadequate approach was used to obtain u(l) and 
further an approximate series expansion was employed f or 
v — 11 [see Eqs. (38) and (39) in lGopakumar fc Iverl J2OO2M . 
As explained earlier, we used a highly accurate and efficient 
method to compute u(l) and an exact relation for v — u, 
given by Eq. 122H . to obtain the time evolution of expres- 
sions appearing in Eqs. 1141 . 

In order to make sure that the numerical codes that 
provided power spectra for H x \q and H+\Q do not con- 
tain bugs, we performed the following consistency check. We 
compared several temporal plots for H x \q(1), created with 
the help of Eq. I|16fl . with those originating from the pro- 
cedure detailed in Section [5] and [3] which does not require 
the evaluation of various Fourier coefficients. For any given 
orbital configuration, we found excellent agreement between 
these two temporal plots, created with two distinct proce- 
dures. A similar analysis was also performed for H+\q(1). 

We note that Eqs. I|16|l and i!2l H . representing Fourier 
series of -ff x |<2(0 and H+\q(1), are the starting points to 
investigate the data analysis issues, associated with h x |q (i) 
and h+\Q(t). This is because the signal received at LISA 
will have induced amplitude, frequency and phase modula- 
tion s on h x \o(t) and h+\ o(t), as detailed in iGutleil (|l998h 
and iRoean fc Bosa i2006T) . It is convenient to incorporate 
these modulations using Eqs. 1161 and 121 H rather than us- 
ing Eqs. Q. This is one of the reasons for explaining in 
detail how we compute Fourier domain versions of Eqs. Q. 

In the next section, we present our results and explain 
their salient features. 



5 DISCUSSIONS 

Let us begin by presenting several plots that depict our com- 
putations. In Figs. 0and|5J we present scaled h x \Q(l) and 
H+\q(1) (H x \q(1) and H+\q(1)) for stellar-mass compact bi- 
naries (mi = rri2 = 1.4 Mq) for various eccentricities when 
n ~ 6.28 x 10 _3 Hz and i = 7r/3. The associated normal- 
ized power spectrum is also displayed in Figs. |3] and 3] We 
clearly see, as expected, as we increase the value of e t , higher 
harmonics with appreciable strengths appear and the total 



power gets distributed among several frequencies. The burst 
nature of the GW signal at high eccentricities indicates the 
strong emission of gravitational waves near the periastron 
due to higher relative velocity around the periastron. An- 
other point to note is that spectral lines are all apparently 
at multiples of radial orbital frequency f r and we do not ob- 
serve the effect of periastron advance in Figs. [3] and 0] Recall 
that the periastron advance causes splitting and shifting of 
spectral lines compared to the Newtonian case. However, in 
our cases, the shifting of lines from their (Newtonian) posi- 
tions is too small to be visible in our graphs. The dimension- 
less parameter k, characterizing the periastron advance, is 
~ 5.9 x 10" 5 for e t = 0.1 and ~3.1x 10~ 4 for e t = 0.9 [For 
n ~ 6.28 x 10" 3 Hz and m = 2.8 M @ ). The frequency shift 
caused by k, deducible from the previous section, is 2k f r and 
in the case of our stellar-mass binaries this is ~ 1.2 x 10 _7 Hz 
for e t = 0.1 and ~ 6.2 x 10~ 7 Hz for e t = 0.9. Therefore, it is 
reasonable to expect that one year LISA observation, which 
can lead to frequency resolution ~ 3 x 10 _8 Hz, should be 
sensitive to k and hence its effect cannot be neglected in 
the search templates. It is also straightforward in our pre- 
scription to drop the effect of k, which can be achieved by 
neglecting 1PN corrections to the orbital dynamics given in 
Section |21 We display Table Q to show that the periastron 
advance creates closely spaced triplets of spectral lines, cor- 
responding to different values of j, which are not visible in 
Figs.|3]and|S] Though triplets are created, one of them really 
dominates the other two in strength. Further, the dominant 
one always corresponds to the lowest possible j value for 
that particular triplet. In another set of plots (Figs. |5] and 
we examine the effect of orbital inclination on H+\q and 
its normalized power spectrum. We observe that the influ- 
ence of i is more visible in the temporal evolution rather 
than in the associated frequency spectrum. We do not dis- 
play similar plots for H x \q as its simple dependence on i, 
defined by cosi, is clearly visible in Eq. Q. We also investi- 
gated if the ratio of the relative power spectrum of H+\q and 
^x|q can be used to obtain information about the orbital 
inclination. From Fig. Q we infer that the above mentioned 
ratio is rather independent of eccentricity and can be used 
to estimate the value of i. 

Let us now discuss why our prescription to compute 
hx\a a nd h+\o is superior to what is available in the lit- 
erature jMoreno-Garrido et al.lll995l . iPierro et alJl200ll) . In 
iMoreno-Garrido et alJ (Il995h . though the orbital motion is 
Newtonian accurate, the effect of periastron advance was in- 
troduced by hand. This was achieved by imposing that the 
periastron advance causes the relative position of the ob- 
server to change in a uniform manner with respect to the 
semi-major axis of the orbit. This leads to certain ad-hoc 
splitting and shifting of the spectral lines. Further, the effect 
of periastron advance and the Newtonian accurate orbital 
dynamics, available in lMoreno-Garrido et alJ dl995l) . do not 
require to specify the individual masses mi and m,2 , but only 
the total mass m. However, our fully 1PN accurate orbital 
description de mands m and 77 an d hence the specification 
of mi and m?. lPierro et alJ (1200 j) . while neglecting the im- 
posed effect that mimics the true periastron advance, em- 
ployed Newtonian accurate orbital motion and an approxi- 
mate analytic expression for the Bessel functions required to 
compute the infinite series appearing in the temporal evo- 
lution for Ii x \q and /i+|q. Recall that J n (ne) and J n (ne), 
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Bessel funct ions of the first kind and the ir firs t derivatives, 
availa ble in iMoreno-Garrido et alJ lll995l) and iPierro et alJ 
(2001), arise from the Fourier analysis of the classical (New- 
tonian accurate) Keplerian motion. 

While comparing our comp u tation s wit h those pre- 
sented in IMoreno-Garrido et al.l (Il995l) and IPierro et al 
20^), we realized that Eqs. fl3~) of | M ^rjng-^^ra^^ ^t^d 



this drift in one year is 



1 199E ) give correct circular limit ( Blanchet 
h x Iq and h+\Q. However, Eqs. (6) and (7) of IPierro et alJ 
(2001) differ, in the ci r cular limit, from Eqs. (2), (3a) and 
(4a) of iBlanchet et alJ dl996fi by a factor of -y/2. The dif- 
ference in sign may be associated with the usage of a differ- 
en t set of conventions c ompared to the one used here and 
in lBlanchet et ail lll996l) . However, this is not desirable be- 
cause for the data analysis purposes, h x \cj and /i+|q need 
to be combined with F x and F+ , the so-called beam-pattern 
functions of the GW detector (Thome 1987J). The conven- 
tion, i.e., the way of prescribing a triad that specifies the 
direction and orienta tion of the binary orbit, is chosen in 
IBlanchet et alJ dl99ri) and here such that the resulting h x \ q 
and H+\q can be dir ectly combined with expressions for F x 
and F+, available in iThornd Jl987h . 

However, as our ap proach is semi-numeri c al, a lmost 
analytic investigati ons of IMoreno-Garrido et al.l il995T) and 
IPierro et alJ (1200 lT) can be used to check our results, es- 
pecially for small eccentriciti es. We restricted comparisons 
of o ur results wi t h tho se of IMoreno-Garrido et all dl995l) 
and IPierro et alJ (l200ll) to moderate eccentricities as it is 
rather difficult and time consuming to compute accurately 
large number of J n (ne) and J n (ne), required for high ec- 
centricities. We find good qualitative agr eeme nt with plots 
availa ble in IMoreno-Garrido et alJ Bl995T> and IPierro et al. 
(2001). Finally, we note that our approach is quite close 
to the way Wahlquist obtained hx\cj(t) and h+\c)(t), with 
Newtonian accurate orbital motion, in the context of space- 
craft Doppler detection of gra vitational waves from widely 
separated massive black holes (IWahlauisdll987l) . It is inter- 
esting to note that the KE is expressed in terms of v rather 
than u in lWahlauisd il987f) and hence inversed trigonomet- 
ric functions appear in the classical KE [It is quite subtle 
to solve the KE expressed in terms of v]. As expected, our 
waveforms with Newtonian accurate or bital motion, are in 
excellent agreement with those given in IWahlauisli lll987D . 

We, therefore, conclude that our fully 1PN accurate 
prescriptions to compute temporal evolutions for Ii x \q and 
/i+|q and their associated power spectra, which do not re- 
quire the infinite series expansions whose coefficients contain 
J nine) and J n (ne), are always going to be more accurate 
and computationally efficient compared to those currently 
available in the literature. 

A rough indicator for the applicability of our waveforms 
can be obtained by comparing LISA's frequency resolution 
for one year, A f lis a ~ 3x 10 _8 Hz, with the frequency shifts 
caused by the advance of periastron and radiation reaction. 
The frequency shift of any harmonic due to the periastron 
advance, as explained earlier, reads 
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The radiation reaction also causes a secular drift in and 



The above expression is obtained by computing the change 
in f T using th e leading order orbit al averaged contributions 
to dn/dt from lDamour et alJ {2004). The rough limits on the 
applicability of our waveforms are provided by investigating 
if LISA can resolve the above two frequency shifts. 

We provide rough limits on the applicability of our 
^x,+ |q(0 by considering three binary configurations. Let us 
first consider the case, where mi = n%2 ~ IAMq, et = 0.9 
and f r ~ 10" 4 Hz. For this binary, Af k ~ 1.3 x 10~ 8 Hz and 
Af RR ~ 4 x 10~ 10 Hz indicating that Newtonian accurate 
orbital description will be sufficient as AfLisA > Af^ and 
Af LISA > Af RR . For NS-NS binaries, with f r ~ 10 _3 Hz, 
e t ~ 0.5, Af k ~ 1.5 x 10- 7 Hz and Af RR ~ 7.5 x 10 _9 Hz 
indicating that our ^x,+ |q(0 with 1PN accurate orbital dy- 
namics will be important to detect these binaries. We do 
not recomm end the use of ad hoc g ravita tional waveforms, 
available in IMoreno-Garrido et al.l dl995f) . as these h x ,+ \Q 
will not provide any information about n and are not com- 
putationally efficient as explained earlier. A rough upper 
limit of applicability is provided by BH-NS binaries with 
m ~ 1OM , e t ~ 0.2 and f T ~ 10~ 3 Hz This is because for 
such binaries 2 x Af RR ~ A f lis a (1 year) and therefore, it 
is important to have /i x ,+|q that include effects of gravita- 
tional radiation reaction. While treating extremely eccentric 
binaries at 1PN order, care should be taken so that does 
not exceed unity. This is also a limit of applicability of our 
h x .+ \g(l) with 1PN accurate orbital motion. 

It should be noted that in all the cases above we as- 
sumed that the emitted gravitational waves, in some cases 
at higher harmonics, are sufficiently strong in amplitudes to 
be observed by LISA. A realistic investigation about the ad- 
vantage and ap plicability of our h x .+ \o(l) is currently under 
investigation in lBose et, al.l <f2C 



6 CONCLUSIONS 

We provided an accurate and computationally efficient pre- 
scription to obtain GW polarizations associated with non- 
spinning bound compact binaries of arbitrary eccentricity 
and mass ratio, moving in slowly precessing orbits. In PN 
description, we restricted the orbital motion to be 1PN 
accurate and considered only quadrupole contributions to 
the amplitudes of h x and h+. We employed an accurate 
(numerical) method due to Mikkola to solve the relevant 
KE appearing in our orbital description. The time domain 
h x \q, h+\g and their associated frequency spectrum, pro- 
vided in this paper, should be required by LISA to detect 
and analyze gravitational waves from stellar-mass detached 
compact binaries moving in eccentric and precessing orbits. 
Therefore, we expect that our accurate and efficient pre- 
scription to compute /i x |q and /i+|q should be of definite 
interest to the recently initiated mock LISA data challenge 
task force. It is possible to extend ou r prescription to in- 
clude effects due to radiation react ion ijDamour et al]l2004l : 
iKonigsdorffer fc Gopakumarl I2006T) ■ spin-orbit interactions 
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( Konigsdo rffer fc GopakumaJ Eooll) and amplitude correc- 
tions. Many of these projects are currently in progress. 
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mean anomaly 1 mean anomaly 1 



Figure 1. Time domain plots of H x \q(1) for various eccen- Figure 2. Plots, similar to Fig. 1, for H+\q(1). 

tricities. The other orbital parameters are mi = m,2 = IAMq, 
% = tt/3 and n = 6.28 X 10 _3 Hz. 
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frequency in units of f 



frequency in units of f 



Figure 3. The normalized relative power spectrum plotted 
against the frequency in units of f r associated with H x \ q dis- 
played in Fig. As eccentricity increases, the dominant har- 
monic shifts its position, and the total power gets distributed 
among the higher 'harmonics'. 



Figure 4. The normalized relative power spectrum associated 
with various H+\q plots displayed in Fig. l2l 
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Figure 5. Plots for H+\q that explore its dependence on the 
orbital inclination. We used the following orbital parameters: 
mi = m 2 = IAMq, et = 0.5 and n ~ 10~ 3 Hz. 



Figure 6. Plots of the normalized power spectrum for various 
values of i associated with Fig. |2J The dependence on i is 
rather difficult in perceive. 
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Figure 7. Plots that depict the ratio of the total power 
present in H x \q and H+\q as functions of the orbital in- 
clination for various eccentricities. The ratio is independent of 
eccentricity and hence can be used to bound i. 
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Table 1. Spectral lines associated with H+\q, and with normalized strengths, appearing at several fre- 
quencies (in units of f r ) for various et- We let mi = ni2 = IAMq, n ~ 6.28 X 10 — 3 Hz, and i = n/3. 
Though the summation index j, appearing in Eq. 1211 . takes integer values, spectral lines appear at \j — 2p\, 
j and |j' + 2p|. We note that different values of j contribute to a closely spaced triplet and one of these lines 
dominates, in its strength, the other two. We observed that normalized value for P$ is always negligible. 
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